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Using heavy meson chiral perturbation theory, we consider the light quark-mass and spatial 
volume dependence of the matrix elements of AB = 2 and AC = 2 four-quark operators relevant 
for B?.-Bf s) and D°-D° mixing, and the B s meson width difference. Our results for these matrix 
elements are obtained in the Nf =2 + 1 partially quenched theory, which becomes full QCD in the 
limit where sea and valence quark masses become equal. They can be used in extrapolation of lattice 
calculations of these matrix elements to the physical light quark masses and to infinite volume. An 
important conclusion of this paper is that the chiral extrapolations for matrix elements of heavy- 
light meson mixing beyond the Standard Model, and those relevant for the B s width difference are 
more complicated than that for the Standard Model mixing matrix elements. 
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I. INTRODUCTION 



Neutral heavy-light meson mixing systems play a crucial role in precision tests of the Standard Model and the 
search for new physics. With the recently measured Am s we can hope to obtain stringent constraints on the 
unitarity triangle of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, provided that the hadronic matrix elements of 
the B°—B° and B^—B^ mixing processes are reliably calculated. On the other hand, the D°—D° mixing system is a 
good channel to search for new physics @], because the Standard Model contribution is strongly suppressed. 

In the Standard Model, the short distance contribution to the mass differences of the heavy neutral meson mixing 
systems (B Q ~B°, B®—B® an d D°—D°) is predominantly determined by the matrix elements of a single set of four 
quark operators: 

O haa = ^(1-75^^(1-75)^, (1) 

where h is a heavy quark field (either a b or a c quark), q a is a light-quark field with flavour a (a is not summed 
over), and a and (3 are colour indices. Models containing flavour-changing currents other than the V — A form (arising 
in supersymmetric extensions of the Standard Model and other scenarios) usually result in mass differences that 
additionally depend on matrix elements of the four-quark operators [|| 
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(the right-handed analogues of Oi^ aa for i = 1,2,3 can also contribute but their matrix elements are the same as those 
above as the strong-interaction conserves parity) . Generically we can represent these operators as 

Oi,aa = hT 1 qhT 2 q 1 (3) 

for the appropriate choice of spin and colour matrices, I\2. In lattice calculations, it is convenient to perform a Fierz 
transformation which renders linear combinations of the operators in Eq. into products of colour-singlet currents. 
We choose to work in the basis of Eq. @ . 

A subset of the operators in Eqs. (QJ and ([2]) are also relevant for calculation of the width-difference in the B s 
system, ATb s /^b b - This difference is the largest amongst the beauty hadrons (Ar£ s /r£ s = 0.31 till [J]) and, 
following an operator product expansion, is given by 

AF m 2 

= I^l^ 7 "' 2 ^ [G(z)(B° s \O hss \B° s ) + Gs(z)(B s \O 2 , ss \B° s )] + 0(1/M b ), (4) 
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where the functions G(z) and Gs{z) are known at NLO in perturbative QCD [f|. At 0(1/Mb), matrix elements of 
C?3 )SS also enter Q- 

Lattice QCD is the only method for calculating the B°—B° and D°—D° matrix elements of the operators in Eqs. ([1]) 
and from first principles, and much effort has gone into such calculations (see Ref. |8| for a recent review). However 
the existing lattice calculations have been performed at light quark masses significantly larger than the physical values 
and necessarily in finite volumes. The effects of these approximations need to be understood. In this paper we consider 
the light-quark mass extrapolation to the physical values for the lattice calculations of these matrix elements. Our 
framework is heavy meson chiral perturbation theory (HMxPT) Ijj j, Hoi . Hlj at finite volume [l2j . The standard model 
AB = 2 operator Ox,aa has been considered in this context in Ref. [12| and here we extend that analysis to the full set 
of operators discussed above. As appropriate for current and foreseeable lattice calculations, we work in the isospin 
limit of SU(3) heavy meson chiral perturbation theory and give results in the SU(6|3) partially-quenched extension. 
Primarily we treat the heavy quark as static throughout but consider the leading effects of the splitting between the 
heavy-light vector and pseudo-scalar mesons. 

An important conclusion of this work is that the chiral extrapolation for matrix elements of 0i, a a is considerably 
less complicated than that for matrix elements of the operators in Eq. ((2|). Generically, the chiral expansion for 
{B\ a) \O x , aa \B\^ takes the form 

{B\ a) \O x>aa \B\ a) ) 7i (1 + L) + analytic terms, (5) 

where 71 is the leading-order low-energy constant (LEC), L denotes the non-analytic one- loop contributions (chiral 
logarithms), and the analytic terms are from the next-to-leading-order counter-terms in the chiral expansion. However, 
for the operators in Eq. ([2]), the chiral expansion has the generic feature: 

(B {s) \O^ aa \B° {s) ) ^ 7i (1+ L) + 7,'L' + analytic terms, (6) 

where i = 2, 3, 4, 5, % and 7 Z ' are unknown leading-order LECs, and L and L' are different one-loop chiral logarithms. 
Again, the analytic terms are from the next-to-leading-order counter-terms in the chiral expansion. The appearance of 
the second non-analytic term complicates the chiral extrapolation in Eq. ^ because an additional unknown parameter 
must be determined. The origin of this complication is discussed in detail in Section ITTT1 

This paper is structured as follows: in Section [II] we briefly discuss heavy meson chiral perturbation theory before 
turning to the inclusion of the four-quark operators in HM%PT in Section Hill We present the results of the next-to- 
lcading order (NLO) light quark mass and lattice volume dependence of the relevant matrix elements in Section IIVI 
before concluding (Section [V} . Various technical details are relegated to the Appendices. 

Whilst this work was being completed, a preprint describing similar work appeared [T3 |. The conclusions of the 
revised version of that work agree with those presented herein, specifically the forms of the chiral extrapolations in 
Eqs. © and ©. 



II. HEAVY MESON CHIRAL PERTURBATION THEORY 

The inclusion of the heavy- light mesons in chiral perturbation theory (HMxPT) was first proposed in Refs. [9I Hoi. [lil] . 
with the generalisation to quenched 1 and partially-quenched theories given in Refs. [HEII. The l/M P (M P is the 
mass of the heavy- li ght pseudo-scalar meson) and chiral corrections were studied by Boyd and Grinstein [TtJ in full 
QCD and by Booth [1J| in quenched QCD. The field appearing in this effective theory is 




where Pa and Pa,}? annihilate pseudo-scalar and vector mesons containing a heavy quark Q and a light anti-quark 
of flavour a. In the heavy particle formalism, such mesons have momentum = A/pv M + fc M with |fc M | <C Mp and 
is the velocity of the particle. Under a heavy quark spin SU(2) transformation S and a generic light-flavour 
transformation U [i.e., U € SU(3) for full QCD and U G SZ7(6|3) for PQQCD (partially-quenched QCD)], 

HP -> SH^UI- (8) 



1 We do not consider the quenched theory here as quenched quantities are unrelated to those in QCD [13 . 
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The conjugate field, which creates heavy-light mesons containing a heavy quark Q and a light anti-quark of flavour 
a, is defined as 

HP = 7°# (Q) S = (j^V + Pi Q)t 7 5 ) (9) 
which transforms under S and U as 

^ Q) — >U ab H { b Q) Sl (10) 

The chiral Lagrangian for the Goldstone particles is 

£gp = ^(s)tr [(0„£t) (fl"E) + E t x + x t S ] , (11) 

where £ = exp(2z$//) is the non-linear Goldstone field, with $ being the matrix containing the standard Goldstone 
fields. We use / = 132 MeV. In this work, we follow the supersymmetric formulation of partially quenched chiral 
perturbation theory [PQxPT] [ljj H3|. Therefore £ transforms linearly under 5t/(3) L <8> St/(3) R and ST/(6|3) L <8> 
£?7(6|3)r, in full QCD and PQQCD respectively. The symbol "(s)tr" in the above equation means "trace" in chiral 
perturbation theory (xPT) and "supertrace" in PQxPT where the flavour group is graded. The variable \ is defined 
as 



X = 2B M q = -^J^I Mqi (i 2) 



-2(0|uu_+ dd\0) 
T 2 

where the quark mass matrix M. q is 

M q QCD) = diag(m u , m n ,m s ), (13) 

in full QCD, and 

A^pqqcd) = diag( ro M , m u , m, , mj,mj,m r , m u ,r n u ,m s ), (14) 

valence sea ghost 

in PQQCD. We keep the strange quark mass different from that of the (degenerate) up and down quarks in the 
valence, sea and ghost sectors. Notice that the flavour singlet state $o = (s)tr(<i>)/v / 6 is rendered heavy by the U(1)a 
anomaly in QCD and PQQCD [3, HH and has been integrated out. 

Furthermore, the Goldstone mesons appear in the HM;\TT Lagrangian via the field 

t = e**'f, (15) 

which transforms as 

£ — Utfrt = U£Ul (16) 

where C/jyr) is an element of the left-handed (right-handed) SU(3) and 5f7(6|3) groups for QCD and PQQCD 
respectively. The HM^PT Lagrangian, to lowest order in the chiral and 1/Mp expansion, for mesons containing a 
heavy quark Q and a light anti-quark of flavour a is then 

£hm x pt = -itTD^W^i^ 

+^tr D (BPhw [?d^-td^] ab HP) 

+5,^7 tr D (ffW)ff«) 7flT5 ) (s)tr - $d^] , (17) 

where B v i = for full QCD, and B„i — 1 for PQQCD 2 . The flavour (super-)trace (s)tr is taken in the appropriate 
flavour space and tro is the trace over Dirac space. The low energy constant (LEC) g occurring in this Lagrangian is 



2 However, since we integrate out the rf in PQQCD 11411 . the coupling 7 does not appear in the results presented in this paper. 
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common to both HM^PT and partially-quenched HMxPT. Note that factors of ^/Mq and ^/Mq have been absorbed 

into the heavy meson fields so the are of mass dimension 3/2. 

The HMxPT Lagrangian for mesons containing a heavy anti-quark Q and a light quark of flavour a is obtained by 
applying the charge conjugation operation to the above Lagrangian [22|]. The field that annihilates such mesons is 

Hp = (P^V-Pi Q h 5 )^, (18) 

which transforms under S and U as 

Hp -^U ab Hl Q) Sl (19) 

The effects of chiral and heavy quark symmetry breakings have been systematically studied at next-to-leading order 
in full [13] and quenched HM^PT [l^]. Amongst them, the only relevant feature necessary for our calculations are 
the shifts to the masses of the heavy-light mesons. These shifts are from the heavy quark spin breaking term 

^tr D (B^HP^) , (20) 

and the chiral symmetry breaking terms 

X 1 B tr D [iM^ + ^M^^H^+X^Bo tr D [iM^ + ^M^]^. (21) 

We choose to use a field redefinition that allows us to work with the effective theory in which the heavy-light pseudo- 
scalar mesons that contain a heavy quark and a u or d valence anti-quark are massless. Notice that the term 
proportional to X[ in Eq. (|2"Tj) causes a universal shift to all the heavy-light meson masses. This means that the 
propagators of the heavy mesons are as follows 

i -iig^u-v^Vu) i ^ -i(9t,u - Vjj ^ 



2(v-k + ie)' 2(v ■ k - A* + ie) ' 2(v ■ k - S us + ie) ' ' 2(v ■ k - A* - S us + ie) 
for P, P* , P S) and P* , respectively. The mass shifts can be written in terms of the couplings in Eqs. (|20|) and ([2lTl : 

A * = - 8 ib> ^ 

and 

S us = 2XiB (m s - m u ). (24) 

In the partially quenched extension, there are two additional mass shifts because the sea quarks masses differ from 
those of the valence and ghost quarks: 

S jr = Mp s - M p = 2XiB (m r - rrij), (25) 

and 

S uj = M p - M P = 2XiB (mj - m u ). (26) 

where P (P s ) is the heavy- light pseudo-scalar meson with a j (r) sea anti-quark. The propagators of the heavy mesons 
containing sea anti-quarks are: 

* -i(a^ - v^Vu) i ^ d —i{j9nv - v^Vu) 



2(v ■ k — d U j + ie) ' 2{v ■ k — A* — 5 u j + ie) ' 2{v ■ k — S U j — 5j r + ie) ' 2(v ■ k — A* — S u j — Sj r + ie) ' 

' ' . (2?) 

for P, P* (vector meson with a j sea anti-quark) , P s , and P* (vector meson with an r sea anti-quark) , respectively. 
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III. FOUR-FERMION OPERATORS IN HEAVY MESON CHIRAL PERTURBATION THEORY 



A. Construction of the AB = 2 and AC = 2 operators 

Under a chiral transformation, the four-quark operators in Eqs. fl} and (|2|) fall into two categories: 

Oll = h T LL q L h T LL q L , 
Olr = h q L h T^n qn, 



(28) 



where 



1±75 



(29) 



Operators Oi tCia , Oi^ aa and 03 jOa are of the first type and transform in the symmetric (6l, 1r) representation built 
from the direct product (3l, 1r) ® (3l, 1r) = (6l, 1r) © (3l, 1r) under chiral rotations while 0±^ aa and 0^^ aa are 
of the second type and transform in the (3l,3r) representation. Here we refer to the SU(3) flavour transformation 
properties, leaving the partially quenched extension to the following subsection. Note that the colour indices in Eq. © 
are relevant to short-distance physics, and hence play no role in the chiral properties of these operators [23] . Treating 

I'll j ^lr an d ^lr as spurions transforming as 



Tll - 


-» S Yll Ul 




P (i) 

1 LR, 


-> <? r (1) 

^ LR, U L' 




r (2) 
1 LR 


^ ° 1 LR u Ri 


(30) 



the operators in Eq. (|28p remain invariant under heavy-quark spin and chiral rotations. We then find that the 
bosonisation of the operators in Eqs. (p} and {2]| is given by 



O 



HM X PT _ 



E 



(fSW)^^] tr D [(zH&>)t-3! x 



(3) + 

aWtr D 



^^r^^^r^] l,(3i) 



for i = 1, 2, 3 where r = Ti = r 2 in Eq. ([3]) and S x and E^. are all possible pairs of Dirac structures 3 . For i = 4, 5 the 
HM^PT operators are 



HM X PT 



E 



-a-^tro 



al 4 2 tr D 



(e^w)^!^] tr D [(Ww) Q r 2 ~; 

(^w^s* (et^) o r 2 s; 
(^w)^!^] tr D [(tHto)r a * 



(^w) a r 2 s x ] tr D [(etffW)^^ 



(6)4. 



tr D 



CO 



(8), 

a^'trn 



(32) 



The positions in the above operators in which the arbitrary Dirac structures, S x and S' x , are inserted is constrained 
by the heavy-quark spin symmetry 0, • Notice that in general both single and double Dirac trace operators must 
be considered. 

Performing the Dirac traces for the particular T, I\2 in Eqs. ([1]) and ([2]), and keeping only the terms that will 
contribute to the matrix elements we consider, leads to the following set of operators involving the individual heavy 



3 An ovcrcomplete list of the possible pairs of structures is: {^x, S^} = {{1, 1}, {1, jj}, {tf, tf}, {7^1,7^}, \siiii>, 7 M }, {_7m^! 7 m ^}> Wfi"''"}, 
{_<?^vi>, \P\ivi>, G^tf}}, their permutations, and possible combinations with 75. There is some redundancy here as the equations 

of motion of the heavy meson fields, ^H^' = etc., relate various terms in Eqs. (13 1 I t and (132 D . 
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meson fields: 



o™* PT = & [(6P (/i)t ) (hp ( v) + (cp; m ) (ep*(^) ' , 

020) {tP m ) a {tP (Jl) ) a +/32(3) {^ m ) a (tP* (n) '») a , (33) 

+^4( 5) (V; (/i)t ) a (ep* ih) '») a +/3i (5 ) (V^O (ep* (S) ^) a , 



(nHMxPT 

2(3), oa 

°4(5), aa 



where the /3-,/3i, and /?■ are linear combinations of the various a~p x appearing in Eqs. l[3"Tj) and ([3"2"| . 

It is important to note that in the above equation, the operator 0^ a * FT behaves somewhat differently from the 
other operators as only a single LEC, occurs. This greatly simplifies any chiral extrapolation of corresponding 
lattice data for neutral heavy-light meson mixing in the Standard Model, as confirmed by the one-loop results presented 
in the next section. We stress that this simplification is not obvious from the operator structure in Eq. [)3ip and is 
particular to the V — A structure of the Standard Model currents. In general, one would expect from Eqs. (|3Tj) and 
([3"2| that operators for pseudo-scalar and vector meson mixing processes are accompanied by different LECs. This is 
the case for all the non-Standard-Model operators, as shown in Eq. ([3"3"[) . 

To understand the origin of the above simplification in the Standard Model operator Of^ xPT , we turn to heavy 
quark effective theory (HQET) [2p, [2?], [Hj]. In this effective theory, the operators that produce the same matrix 
elements as those in Eq.@ are [29[ 

^hqet nr-n nir„n -i-nt 



i. a a 



QT iqa Q^T 2 q a + Q^T iqa QT 2 q ai (34) 



where F1.2 are the appropriate Dirac and colour structures from Eq. (T5|). Here, Q and Q denote fields annihilating a 
heavy quark and heavy anti-quark, respectively (these fields do not create the corresponding anti-particles). Additional 
terms in HQET which create two heavy quarks or annihilate two heavy anti-quarks will not contribute to neutral 
heavy-meson mixing and are ignored. 

The standard model operator in HQET, Of „<f T , satisfies the relation 

K«r T } i-p) = {ss.ogr} 1^ = 0, (35) 

where \P) is pseudo-scalar heavy-light meson state, and \P) is the state of its anti-particle. The operator 

Sq = <' ,:i (y*,,(] - Q<r iS &], (36) 

is the heavy quark spin operator [30] that changes the spin of the heavy-light meson state by one. Therefore, Eq. ([331) 
implies that the mixing matrix elements for vector and pseudo-scalar heavy-light mesons are equal and opposite [221 ] 
in the heavy-quark limit. This symmetry is reflected in HM%PT, leading to the result for T in Eq. ([33[) . 

For the non-Standard-Model operators, it is straightforward to show that 

{S 3 Q , « ET } \ p ) * 0, {s%, O^) \P) * 0, (37) 



and 



o3 //-)HQET 



P)+% SL0f^ T \P)?0, (38) 



where i = 2,3, 4, 5. This means that the pseudo-scalar and vector meson mixing processes via these operators are not 
proportional to each other, hence the appearance of the terms accompanied by [3' 2 345 and (3± 5 in Eq. (j33[) . 

We end this subsection by noting that equations of motion for the heavy quark [3l| result in O^^a = —Oi >aa /2 — 
2 aa, and can be used to relate some of the LECs in Eq. 



B. Partially-quenched extensions 



In the partially-quenched theory, the operator matching is analogous because the QCD operators considered on 
the lattice involve only valence quarks. Since HM^PT is contained within partially-quenched HM^PT, the LECs 
occurring in the four-quark operators of both theories are the same for the quantities we consider. 
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(a) (b) (c) 

FIG. 1: Diagrams contributing to the matrix elements of four-quark operators at NLO in the chiral expansion. Solid, double 
and dashed lines correspond to propagators of pseudo-scalar and vector heavy-light mesons, and Goldstone mesons, respectively. 
The crossed circle denotes the four-quark operator and diagram (a) is the wave-function renormalisation. 



In the partially quenched case, the AB — 2 and AC = 2 operators transform in the symmetric tensor product of 
two fundamental representations of SU(6|3), a 42-dimensional representation. The operators arising from QCD can 
be simply embedded in this larger representation with no mixing into different representations. In most cases (and 
herein) it is sensible to choose the quark "charges" such that the operators are purely valence, but any other clement 
of this representation suffices. 



IV. NEUTRAL MESON MIXING MATRIX ELEMENTS 



Calculations at NLO in the chiral expansion require the evaluation of the one-loop diagrams shown in Figure [TJ 
We perform these calculations both at infinite volume and in a cubic spatial box of dimensions L 3 (the time extent is 
assumed to be infinite). In this section we summarise the results, relegating details of the calculations to Appendix [Al 

For the standard model operator we find the following matrix elements 
(B°\O hdd \B°) = (3, 



(B° s \O hss \B° s ) = (3, 

The wave-function contributions, Wat, and tadpole- and sunset-type operator renormalisations, Ta an d Qa^ (dia- 
grams (b) and (c) in Figure [TJ respectively) are non-analytic functions of the light quark mass and lattice volume 
and are defined in Appendix |5J The "analytic terms" here include Goldstone meson mass squared terms, a term 
~ a s (Mf,) /47T (arising from mixing) and a term ~ Aqcd/M^. The ~ Aqcd/-W& term arises from l/Mj, terms in the 
Lagrangian (those from Eq. (f2T)]) are included) and from additional 1 /M^-suppressed HM vPT operators that match 
onto the QCD operators. We note that at higher orders, the simple relation between the B and B* matrix elements of 
the Standard Model operator will break down. Although we present results specifically in the B-meson systems, note 
that they are also applicable to D-meson systems, under the assumption that the charm-quark mass is large enough 
compared to Aqcd- 

Parameterising the AB = 2 matrix elements in the standard form [3[ (/b is the weak decay constant defined 
through (0\hnJ5q\Bq(p)) = ip»fj3 q ), 

{B Q q \0 lm \B») = \mIjIB^), (40) 
(B° q \O im \B° q ) = r H R 2 M 2 Bq f 2 Bq B^) for i = 2, . . . , 5 , (41) 

(R = — fJ^ra TTIT an< ^ ^ 2 = ~ I' ^ 3 = i> ^ 4 = ^' ^ 5 = f) ^ uc parameters, B B \/s) agree with those derived in 
Refs. [fill where the relevant expressions for fg are also provided. 

For the additional operators that contribute to the S-meson mixing processes beyond the Standard Model, we 



1 + T« + Ws " ± Wb " + ) + analytic terms , 



(39) 



1 + T} 1 ' + "- + QW + analytic terms . 
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obtain: 



(B \O 2(3Ud \B°) 



(B° s \O m , ss \B° s ) = 



(-B°|04(5),dd|-B°) 



$ 2 (3) (l + rr )] + W * 0+ 2 Wb ^ + ^ (3) 3 » + analytic terms , 
$2(3) i 1 + T i 2(3)) + Wb " + 2 Wb °^ + /3 2(3) ?<»» + analytic terms . 



(B° s \O M5) , ss \B s ) = 



$4(5) + $4(5) 
$4(5) + $4(5) 



(4(5)) Wgo + VVgo 



i + r: 



1 + T^ 5 " + 



)+ 


$4(5) " 


1" $4(5) 


)+ 


$4(5) " 


1" $4(5) 



(42) 



Q d + analytic terms , 



The terms ~ Q q arising from the sunset diagrams [Fig. Hfc)] involve the neutral heavy-light vector meson mixing 
amplitudes. As discussed in the preceding section, it is only in the case of 0\ m that these amplitudes are related to 
those of the pseudo-scalar heavy-light mesons. For i = 2, 3, 4, 5 these terms are consequently accompanied by different 
LECs. The analytic terms in the above expressions depend on the renormalisation scale is such a way as to cancel 
the scale dependence of the non-analytic loop contributions. 

The matrix elements in Eqs. ([59]) and (|4"2"| also determine the B s decay width differences [Eq Q]. The extrapolation 
in light quark mass and lattice volume is more involved here than for the matrix elements determining Standard Model 
oscillations and have not been accounted for in the existing unquenched lattice calculations [32j,[33j. Direct calculations 
of the ratio of the matrix element of Oi jSS to that of 02, ss do not help in this regard as the non-analytic behaviour 
does not simplify. 

At present, we can only study the finite volume and light quark mass effects described in these formulae cursorily 
as there is very little lattice data to use for such a task in any reliable manner. This will hopefully change in the near 
future; the recent calculations of (33 | are encouraging (we note, however, that our results imply that more than three 
light quark masses are needed for the NLO chiral extrapolation of matrix elements of 0i, a a for i = 2,3,4, 5). As a 
guide to the importance of such effects, we present representative results for the 02,4 matrix elements. Results for the 
Standard Model operators have been discussed in Ref. [13]. Taking the QCD limit for definiteness, Figures and [3] 
explore the dependence of the finite volume shifts in the matrix elements of 02, dd, 02, ss 04, dd and 04, ss , normalised 
by their tree-level values, on the pion mass for two different volumes, L — 2.5, 3.5 fm. In each figure we plot the ratio 



(B° {s} \O l j f \B l 



\ 

(s)/FV 



af {s) \O lJf \B° s) )(L) - (Sj. ) |O il/ /|^, ) )(oo) 



(43) 



(s)' 



We fix / = 0.132 GeV and A» = 50 MeV (variation with A* is small and similar to that found for the Standard 
Model operator pjj], with the FV effect decreasing with increasing A*). Using recent CLEO measurements [35l |36| . 
the coupling g is taken to be 0.3 < g 2 < 0.5 with a central value of g 2 = 0.4; variation in g is indicated in the figures 
by the inner (darker) shaded regions. For each quantity, we vary the ratio of B* to B LECs over a reasonable range, 
taking |$ 2 /$ 2 | < 2 and \{(3' A + $4)/ ($4 + $4)! < 2 (naturalness would suggest these ratios should be of order unity, 
but for simplicity we also allow smaller values). In each figure, the central curve corresponds to a ratio of unity and 
the outer (lighter) shaded region to this variation. As can be seen, effects of the finite volume on Oi.dd are similar in 
size to those found for the Standard Model operator in Ref. pl| • Effects for the strange operators are considerably 
suppressed as pion loops do not contribute to these matrix elements in the QCD limit. 

Before concluding, it is useful to consider how lattice spacing artefacts will enter the above expressions. Since 
calculations here involve both light and heavy quarks, there are two types of discretisation effects. The effects of 
light quark discretisation are very simple to incorporate; at a particular lattice spacing, a, the masses of the various 
Goldstone mesons in chiral loops are shifted from their continuum values (and these shifted masses should be used in in 
fits to lattice data using the above partially-quenched expressions) and the various counter-terms become polynomials 
in the lattice spacing. In general this polynomial will contain all powers of a, but if both the light-quark action and 
the four quark operator are improved (or a discretisation satisfying the Ginsparg- Wilson relation [33] is used for the 
valence quarks (38|), the leading corrections in a can be eliminated. For most foreseeable calculations, this is the 
extent of discretisation effects at the order to which we have worked (we assume that a Aqcd ^5 w 9 /Aqcd)- However, 
if a heavy quark action is used that breaks heavy quark spin symmetry at 0(a), additional complications will arise as 
this symmetry can no longer be used to constrain the form of the EFT operator, 0f^ xPT ; B^-B^ and B* s ^-B^ 
matrix elements are no longer related. The resulting mass and volume dependence of matrix elements of this operator 
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FIG. 2: Finite volume effects in mixing matrix elements of the operators Oi.dd for i = 2,4 for two lattice volumes, L = 2.5, 3.5 fm. 
The central curve correpsonds to g — 0.4 and f3'i/{3i = 1 while the inner (darker) and outer (lighter) shaded regions correspond 
to variation of 0.3 < g < 0.5 and \j3[/j3i\ < 2. The curves terminate at L = 2.5 where p-regime chiral perturbation theory 
becomes unreliable. 
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FIG. 3: Finite volume effects in mixing matrix elements of the operators Oi t3S for i — 2,4 for two lattice volumes, L — 2.5, 3.5 fm. 
Details are as in Fig. [5] 

will become more complicated, resembling instead that of the matrix elements of the other operators in Eqs. pTj) and 

(EH)- 4 



CONCLUSION 



We have considered the matrix elements of four-quark operators relevant for heavy-light neutral meson mixing 
and decay width differences in partially-quenched, finite volume heavy meson chiral perturbation theory relevant for 
lattice computations. For heavy-light neutral meson mixing, inclusion of operators beyond those in the Standard Model 
complicates the chiral extrapolation as two LECs appear at leading-order rather than one for the Standard Model 
operator. The matrix elements relevant for lifetime ratios and decay width differences have similarly complicated 
light-quark mass and spatial volume dependencies. This arises due to the fact the operators for these processes do not 
guarantee that B—B and B*—B* mixing amplitudes are proportional to each other in the heavy-quark limit. Our 
results are useful for current and future lattice calculations of these matrix elements, which are needed in high-precision 
tests of the Standard Model and the search for new physics. 



4 Additional complications beyond the scope of this work may arise if the Kogut-Susskind fermion action is used for the light quarks. 



10 



Acknowledgments 

The authors thank M. J. Savage, S. R. Sharpe and M. Wingate for useful discussions. This work is supported by 
the U.S. Department of Energy via grant DE-FG03-97ER41014. CJDL is grateful for the hospitality of DAMTP, 
University of Cambridge and CTS, National Taiwan University. 



APPENDIX A: INTEGRALS AND SUMS 



We have regularised ultra-violet divergences that appear in the various loop integrals using dimensional regularisa- 
tion, and subtracted the term 



A 



4-d 



7_e +log(47r) + 1. 



The integrals appearing in the full QCD calculation are defined by 

d d k 1 



I- X {m) = n 



4-d 



(2ir) d k 2 -m 2 + ie 



im 

16^2 



A — loe 



(Al) 



(A2) 



H- x (m,A) = («r 

d f d d k 



dAJ {2ir) d (k 2 -m 2 + ie){vk- A + ie) 
d 



where 



F- x (m,A) = 



2A 2 



. /10A 2 4m 2 \ . 2(A 2 -m 2 ) /A 



with 



R(x) = \/ x 2 — 1 log 



^/x 2 — 1 + ie 



x + v 'x 2 — 1 + ie J 

and /z is the renormalisation scale. For the partially quenched calculations, we also need the integrals 

d d k 1 _ dl- x (m) 



{2v) d {k 2 -m 2 + ief 



dm 2 ' 



and 



fl?'(m,A) ee (gP" - »'/) M 4 - d 



d d k 



_d_ 

OA J (2ir) d (k 2 -m 2 + ie) 2 (vk-A + ie) 



() 



dm 2 



H- x (m,A). 



(A3) 



(A4) 



(A5) 



(A6) 



In a cubic spatial box of side length L with periodic boundary condition, the three-momenta are quantised as 

and one instead obtains the sums (after subtracting the ultra-violet divergences) 

J ("0 55 h X- 1 ^ fc 2 - m 2 + ie = 1{m) + JFv(m) ' 



(A7) 



(A8) 



(A9) 
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and 



for the full QCD calculation, where 



(k 2 — to 2 + ie)(v ■ k — A + ie) 



Hm) = h( m )\\=o> 



H(m,A) + H FY (m,A) (A10) 



(All) 



and 



H(m)=H x (m,A)\ x=0 , 
are the infinite volume limits of X and Ti, and (n = \n\) 



(A12) 



hv(m) 



4tt 2 ^ V 2nL V nL 



-nmL 



X <M + 



15 



8titoL 128(nmi) 2 







nmL 



(A13) 



is the finite volume correction to /(to). The function i?FV is the finite volume correction to H(m, A) and can be 
obtained via 



flpvKA) = 



(to 2 - A 2 )K FV (m, A) - 2A J FV (m, A) + U FY (m) 



(A14) 



where J F y(m, A) and _?Cfv(w, A) are given by 



wr = y |fc| 2 + to 5 



J PV (m,A)= (-) £ / d\k\ 



sin(|fc||n|£) 



and 



if FV (m,A) 



w fe + A ) / V 



dJ F y(m, A) 
dA 



In the asymptotic limit where mL 1 it can be shown that (with n=\n 

JfvKA) = W-i 



n^O 



8tttiL 



e~ nmL A, 



(A15) 



(A16) 



(A17) 



where 



-4 = e^ 2) [1 - Erf (z) ] 



/ 1 



1/2 z 6 \ z A ,,2 



7T 14 2 



nmL 



\nmL 

1 /9z 5z 3 7z 5 z 7 



-e( z ) [l - Erf (z) ] 



y + y)e^[l-Erf(z)] 





1 








1 


nmL 





with 



nmL 



, (A18) 



(A19) 



The quantity A is the alteration of finite volume effects due to the presence of a non-zero A. See Ref. [12| for further 
discussion. 
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For the PQxPT calculations, one also needs 



1 



dl(m) dlFy(m) 



(k 2 — m 2 + ie) 2 



dm 2 



dm 2 



(A20) 



and 



W (m, A) 



<9A 



dH(m,A) dH FV {m,A) 



dm 2 



dm 2 



(k 2 — m 2 + ie) 2 (v ■ k — A + ie) 



(A21) 



APPENDIX B: LOOP CONTRIBUTIONS AND BAG PARAMETERS 



In this appendix, we present results for the various contributions in Eqs. (j39|) and (|42|) . Wb° , Wgo , and 



d(s) 



qW s) (j = 1,2,3,4,5). These results are given in the sea and valence isospin limit of SU(6|3) partially-quenched 
HM^PT, with the quark masses given in Eq. (|14[) . The QCD limit, where sea and valence quark masses are equal, is 
easily taken by setting m,j — m u and m r — m s . We also present the bag parameters defined in Eqs. (|40[) and (|41[) . 



1. Loop contributions in SU(6 3) partially-quenched heavy meson chiral perturbation theory 

To compactly express the partially quenched expressions, it is useful to define the following quantities: 



2(6 2 VS -M 2 + M 2 X ) & 



vs 



(M*-M%) 2 



(Bl) 



R — 1 — A 



(B2) 



r 3/f 2 - 2S ^ S 



(B3) 



A. 



(s^ 



VSs 



{26 2 vs - Ml + M? s ) 



(25 2 vs + AS 2 VSs -M 2 + Ml s y 



^s,s 1 ^S,S 1 



(B4) 
(B5) 



Cs.s — 



2 ) 

s.s ) 



(26 2 vs + A6 2 VSs -M 2 + Ml l 



(B6) 



with M l b = M^a + m b ), 5 2 vs = Ml - M*j, S 2 /Ss 



Ml 



M? r , M„ = M UtU , and M% = f (M» + 2M 2 S 



25 



vs 



For ease of use, we note that in the QCD limit (setting valence and sea masses to be identical), 



A QCD = £ 

u,u 2 ' 



B QCD 



r QCD _ n 
L'„ ,, — u, 



aQCD _ o 



B 



QCD _ 



= -2, 



C?? D = 0. 



The various loop contributions can then be written as 



■ 2 



2 

w B o = w B o = -i^j \b u , u h (m x , a*) - en {m uJ , a* + ^ ) - 3^ (m u , p , a, + <5 sr + 5 US ) 

(M„,„, A*) + C„,„H'' (M» >u , A 
B S:S H (Mx, A*) - (Af SJ -, A, + - M - 3H (M,, r , A, + <5 sr ) 

+A S , S H (M s , s , A,) + C^W"' (M a ,„ A.) 

for the wave- function renormalisations, 

2B U , U 1{M X ) - 6T(M uJ ) - 3T(M u , r ) + (2A M ,„ - 3)T(M U ,„) + 2C U ,J^' (M„,„) 
2B S>S X(M X ) - 6T(M a j) - 3T(M s , r ) + (2A S , S - 3)Z(M s , a ) + 2C. .7"' (M S)S ) 



T (l,2,3) = _J_ 

T (l,2,3) = _J_ 
3/ 2 



^) = -^[2X(M„, i)+ X(M„ )r )-T(M„,„)] 
^ s (4 ' 5) = — ^[2T(M a j)+T(Af., r )-T(M,,.)] , 



for tadpole integrals [Fig Hfb)] and, 

■ 2 

Q( d = Jp ( B u,uH (Mx, A) + (A M , U - 3) H [M u<u , A) + C u , tl W' [M u ,„, A] 

■ 2 

= ^(b s , s H(Mx,A) + (A s , s -3)W(M s , s ,A) + C s H , ''[M s , s ,A] 
for "sunset" integrals [Fig [He)]. 

2. Bag parameters 

For completeness, the bag parameters defined in Eqs. ((40]) and (f4Tj) are given by: 

\ _ 3/3i / g 2 X H , u 



siJjW = # (i - — - — 

(2/3), x fti/3 ( X I>U g 2 (3' 2/3 X H .u 

R (2/3)/,x _ / ,g 2 /? 2/3 ^H, s 

R (4/5), v _ &/5+&/B (, , ^ 5 2 (^/5+^/ 5 )^ 



/34/5+A/5 A , X Z)S , g 2 (/%/ 5 + /%/ 5 ) X H ,. 



2 ^ R2 \ f2 / 2 (/34/5+/?4/5) 
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where k is the LEC governing the heavy-light axial current [12| and for convenience we have defined 




(B21) 



(B22) 



(B23) 



(B24) 
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